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Abstract 

We perform one- and two-points magnitude cumulant analysis of one-dimensional longitudinal 
velocity profiles stemming from three different experimental set-ups and covering a broad range of 
Taylor scaled Reynolds numbers from R\ = 89 to 2500. While the first-order cumulant behavior 
is found to strongly depend on Reynolds number and experimental conditions, the second-order 
cumulant and the magnitude connected correlation functions are shown to display respectively 
universal scale and space-lag behavior. Despite the fact that the Extended Self-Similarity (ESS) 
hypothesis is not consistent with these findings, when extrapolating our results to the limit of 
infinite Reynolds number, one confirms the validity of the log-normal multifractal description of 
the intermittency phenomenon with a well defined intermittency parameter C2 = 0.025 ± 0.003. 
But the convergence to zero of the magnitude connected correlation functions casts doubt on the 
asymptotic existence of an underlying multiplicative cascading spatial structure. 

PACS numbers: 47.27.Eq, 02.50.-r, 47.27.Jv, 47.53. +11 
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Since Richardson's original work [ij , a common "mental image" of fully developed turbu- 
lence is a dynamical cascading process in which large eddies split up into smaller ones which 
themselves blow up into even smaller ones and so forth 0. According to this picture, en- 
ergy propagates from the integral scale, where eddies are generated, down to the dissipative 
scale, where they vanish by viscous dissipation, through a multiplicative process, each eddy 
inheriting a fraction of its parent's energy. Since this early intuitive description, the notion 
of cascade has remained the creed of many models proposed in the literature || to mimic 
the statistical properties of turbulent signals. Among the available experimental data, a lot 
of effort has been devoted to the study of the longitudinal velocity component recorded in 
directional flows such as jets or wind tunnels ||, [| [|, ||. In these configurations, the Taylor 
hypothesis that considers the spatial structure of the flow as globally advected upon the 
probe, enables us to interpret temporal time series as spatial profiles. In 1941, Kolmogorov 
resumed Richardson's picture in his statistical analysis of the spatial fluctuations of ve- 
locity profile in the sense that he linked the one-point statistics of the velocity increments 
5vi = v(x + l) — v(x) over different distances I, by some dimensional analysis which predicted 
the remarkable scaling behavior of the moments of Svi : M q j =< 5v? >~ l^ q where Q — <?/3. 
Actually, ( q turned out to be a non-linear function of q in most experiments [0, ^, f|, [| and 
many studies inspired from Kolmogorov and Obukhov second theory [[/J tried to explain 
and to predict the analytical shape of this non-linearity. The controversial situation at the 
origin of some disagreement between models (e.g. the log-normal [171 |SL M as opposed to 
log-Poisson models [pTo| ) results from the experimental observation that the moments M q j 



do not really scale perfectly. Indeed, there is a persistent curvature when one plots ln(M q j) 
vs ln(7), which means that, rigorously speaking, there is no scale invariance. In order to give 
a sense to the exponents ( q , Benzi et al. |ll| defined the "Extended Self-Similarity" (ESS) 
hypothesis by proposing the following behavior for the moments of the velocity increments : 

M q ,i ~ f(l) c «, (1) 

where /(/) would be some g-independent function of I. Along this line, ln(M gi ;) vs ln(M ?2; /) 
curves look definitely more linear than using standard log-scale representations and, by 
assuming that ^3 = 1 0, some experimental consensus has apparently been reached on the 
nonlinearity of the ( q spectrum . In a recent theoretical work |12| , Arad et al. suggest that 
at low Reynolds number, structure functions scaling properties are "poluted" by anisotropic 
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effects that can be mastered using the irreductible representations of the rotation group. 
Unfortunately, this analysis is not tractable with single point data. Hopefully, as shown 
below, some statistical quantities turn out to display universal behavior that are likely to 
be insensitive to anisotropic effects. 

In the early nineties, Castaing et al. || recasted the cascade picture and the ESS hy- 
pothesis in a probabilistic description that accounts for the continuous deformation of the 
probability density functions (pdf) of 5vi with I, by the mean of a propagator Gw (/' > I): 

/oo 
G w (u)e- u P v (e- u 5v)du. (2) 
-oo 

This equation can be related to a cascade process in which the variable 5vi is continuously 
decomposed as Svi = Yl7=i ^k+iM^i'i where lndW^^/J) are independent random variables 
of law Gi i+1 i r Within this framework, the ( q and the /(/) ESS functions can be related to 
the shape of Gw thanks to the cumulant generating function of G [|, ||, 

hiGwi-iq) = In {M q jM q , v ) = ( q (/(/) - /(/')) , (3) 

the classical scale invariant case corresponding to /(/) = ln(7). In this paradigm, the suc- 
cessive terms of the polynomial development of In (M q ^/ M q ^i) as a function of q, involve 
the cumulants C q (l) = —C q f(l) of the logarithm of \5vi\ ("magnitude") from which one can 
express the ( q spectrum as ( q = — J2t=i ^k%- 

The question whether turbulent velocity signals do (or do not) present a space-scale 
cascade structure thus amounts to checking that the successive velocity magnitude cumu- 
lants possess the same scale behavior. This perspective can be extended from structure 
functions to unfused correlations functions |TJ| (e.g., < 5vf(x)5vf,(x + Ax) >), by study- 



ing the behavior of the magnitude connected correlation functions (MCCF) : Cw{Ax) = 
< \n(\Svi(x)\) ln(|5f/'(x + Ax)|) > c . These correlation functions have been introduced in 
Ref. |1| and play the same role for the multivariate multi-point velocity law as the previous 
magnitude cumulants for the one-point law. For multiplicative cascade processes, they are 
expected to display a logarithmic dependence on the spatial distance Ax ||15|| . Our goal in 
this Letter is to push further the analysis carried out in Refs || by studying both the 
cumulants and the connected correlation functions associated to ln(|<5t>;|) at different scales 
/, with the aim to establish a clear diagnostic about (i) the scaling properties of experimental 
turbulent velocity records and (ii) the validity of multiplicative cascade models in the limit 



of infinite Reynolds number. We will work with signals gracefully supplied by three different 
experimental groups. The signals labelled with Taylor scaled Reynolds number R\ = 89, 
208, 463, 703 and 929 stem from Castaing's group and were recorded in a gaseous helium jet 
at very low temperature in CRTB (Grenoble) |TBJ; the signal labelled R\ = 570 stems from 
Baudet's group and was recorded in an air jet in LEGI (Grenoble); the signal labelled R\ = 
2500 was recorded at the ONERA wind tunnel in Modane by Gagne and collaborators |T7 



Let us start investigating one-point statistics via the computation of the three first cu- 
mulants C q (l) of ln(|<5t>;|). In Fig. la, C*i(Z) is plotted vs ln(Z) for all signals. They are all 
characterized by an integral scale L corresponding to a decorrelation of the increments | Svi | 
for Z > L. Below this saturated regime, there is no well defined "inertial range" since one 
observes a continuous cross-over across scales towards the smooth dissipative regime down 
the Kolmogorov scale 77. Let us note that the integral scale does not depend on the Reynolds 
number but only on geometrical considerations. As an illustration, one gets the same value 
for L for all the signals stemming from Castaing's group which were recorded in the same 
apparatus, R\ being controlled either by the mean speed of the flow or by the value of viscos- 
ity that is changed by tuning temperature. Furthermore, the numerical values obtained for 
L : 1.82 cm for Castaing's group, 0.91 m for Baudet's group and 10.2 m for Gagne's group, 
correspond to the characteristic sizes of the experimental set-ups. L is thus a geometric 
parameter and when one increases R\, the inertial range increases via the decrease of 77. 
Let us notice that as proposed in Ref. JRJ, the shape of the cross-over from the saturated 
regime to the inertial one is rather well reproduced by a simple Langevin equation (inset of 
Fig. Id) : 5v(x) = —jv(x) + aW(x), where 7 = 1/L and W(x) is an appropriate fractional 
Brownian noise (Bh = 1/3). However the precise functional dependence of C\(l) is likely to 
depend on the experimental parameters and to be strongly influenced by anistropic effects 

E3- 

The situation is very different for the second magnitude cumulant C^O when plotted vs 
ln(Z) as shown in Fig. lb. We first notice that the decorrelation scale where C^Z) saturates, 
is even better defined with C^Z) than with C\(l), since a linear behavior is observed up 
to this scale without any cross-over. The most important feature is that, whatever the 
set-up and whatever R\, all the data fall on linear curves which all have the same slope 
— C2 = —0.025 ± 0.003. As a matter of fact, if we superimpose all the curves by setting all 
L to 1, as shown in Fig. le, we cannot distinguish anymore one curve from the others. This 
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consideration is particularly relevant in the case of the signals of Castaing's group where R\ 
is varied without changing L. 

The results concerning the third magnitude cumulant C%(1) are reported in Figs lc and 
If. The slope of the curves obtained when plotting C^{1) vs ln(Z) systematically decreases 
when increasing R\. For small R\, C3 is significantly different from zero which means that 
the log- normal paradigm is not valid. In fact, from the curvature of the corresponding 
experimental C\(l) curves, we believe that this is rather a confirmation of the absence of a 
well defined inertial range for these low values of R\. For the largest values of R\ (> 800), 
C3 becomes small enough (up to finite sample effects as previously reported in Ref. |§) 



to be neglected ]T8[. At high Reynolds numbers, one can thus suppose that C 3 = (and 
thus so the higher order cumulants) that implies a normal shape for the propagator Gw in 
agreement with log-normal models |7], |^, |§. To summarize, the scale behavior of structure 
functions is well described by the first two terms in the cumulant expansion of Gw and 
the differences observed in the behavior of C*i(/) (Figs la and Id) and C 2 (l) (Figs lb and 
le) bring the experimental demonstration of the inconsistency of the ESS hypothesis which 
rigorously requires an identical scale behavior for the two cumulants (Eq. (|3|)). The full 
Cj spectrum cannot thus be defined in the range of Reynolds numbers we have investigated 
but the scaling exponent C2 = 0.025 ± 0.003 of is likely to be a universal (model 

independent) characteristic of intermittency. 

Let us now focus on the MCCF : C u >(Ax) =< \n(\5vi(x)\) \n(\5vi>(x + Ax)|) > c . As 
pointed out in Ref. |15| , we have checked that the MCCF computed at different scales / 
and I' all collapse on a single curve provided Ax > max(/,/'). Let us recall that for a 
multiplicative cascade, we expect the MCCF to behave as — C2 ln(Ax/L). In Fig. 2 are shown 
the MCCF computed at a single small scale (I = /') for each experimental signal. Three 
main observations must be raised from these curves. First, the MCCF do not behave linearly 
but as the square of the logarithm of the spatial distance Ax : Cu(Ax) = a(R\) ln 2 (Ax/L). 
This quadratic behavior is universal as far as experimental set-ups are concerned and has 
never been observed before. Second, the respective integral scales at which the magnitudes 
are decorrelated are nearly the same as the ones observed on the corresponding cumulants 
C q (l) (q — 1;2), which is not a trivial result. Finally, as shown in Fig. 3, the prefactor a(R\) 
has a systematic decreasing behavior as a function of R\. As far as the analytic shape 
of this decrease is concerned, there is no certainty. As illustrated in Fig. 3a, a power-law 
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behavior : a ~ R^ a with a ~ 0.20 ± 0.03 provides a reasonable good fit of the data. But as 
shown in Fig. 3b, one can also fit the data by a ~ K/ ln(R\) with K = 0.027 ±0.006, which 
is a particular Reynolds number dependance that allows us to define a characteristic scale 
of the flow under study. Indeed, if one considers that the decorrelation scale is the same 
for the magnitude cumulants C q (l) and the MCCF, from the behavior of the variance as 
— C*2 ln(7/L) and of the MCCF as aln 2 (Ax/L), one can define a characteristic scale l c where 
the two curves meet (when identifying scale and space-lag according to the multiplicative 
cascade picture |15| ) : ln(/ c /L) = —C 2 /a. With the numerical values extracted from Figs 
lb, 2 and 3b, one gets (l c / L) ~ R^ b with b = 0.91 ± 0.3. This new scale l c might well be the 
Taylor scale A for which 6 = 1. In that respect, our approach could be a way to measure 
objectively A and this even when the dissipative scale r\ is not resolved. This universal 
quadratic behavior of the MCCF vs ln(Ax) and the systematic decrease of a(R\) when 
increasing R\, are new observations which enables us to make some conjecture about the 



asymptotic limit of infinite Reynolds number. As pointed out in Refs ||, [15], correlations in 
the magnitude are the signature of a multiplicative structure. If we extrapolate the observed 
behavior to infinite R\, we predict the convergence to zero of the MCCF on any finite range 
of scales. That is to say, in the limit of infinite R\, the intermittent spatial fluctuations of 
the longitudinal velocity component are not likely to display a multiplicative hierarchical 
structure. 

To summarize, we have advocated in this paper the study of scaling properties of longi- 
tudinal velocity increments by means of one- and two-points magnitude cumulant analysis. 
The results of our measurements for seven flows, at seven different R\ and stemming from 
different experimental set-ups, are two-fold. Concerning the one-point (fused) statistics, we 
mainly observe that the log-normal approximation is pertinent at sufficiently high Reynolds 
number and that the first two cumulants have not the same scale behavior as assumed by 
ESS hypothesis. While C\{1) displays some Reynolds and set-up dependent departure from 
scale invariance, C2(/) exhibits scale invariance behavior with a universal intermittency co- 
efficient C*2 = 0.025 ± 0.003. Concerning the two-points (unfused) statistics, we observe 
that the behavior of the MCCF vs ln(Ax) is quadratic whatever the considered experiment 
but with a prefactor a(R\) which decreases when increasing R\. These new observations 
lead us to conjecture that in the limit of infinite Reynolds number, the space-scale structure 
underlying the fluctuations of the longitudinal velocity is not of a multiplicative type but 
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more likely in the spirit of the log-normal multifractal description pioneered by Kolmogorov 
and Obukhov in 1962 |7|, i.e. with an intermittency coefficient C2 = 0.025 ± 0.003 but 
without any correlations across scales. This is not at all shocking from a physical point of 
view since the "mental image" of a dynamical cascading process does not a priori imply 
that there should be some multiplicative (cascading) spatial organization. Furthermore, one 
must realize that our conclusions come out from the study of ID cuts of the velocity field 
only. In a forthcoming publication, we plan to carry out a similar analysis of a 3D velocity 
field issued from direct numerical simulations, with the specific goal of testing the validity 
of this non multiplicative log-normal multifractal picture to account for the intermittent 
nature of fully developed turbulent 3D velocity fields. 

We are very grateful to B.Castaing's, C.Baudet's and Y.Gagne's groups for the permission 
to use their experimental signals. 
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FIG. 1: Magnitude cumulants C q (l) of the seven studied signals for R\ = 89 (o), 208 (•), 463 (A), 
570 (x), 703 (A), 929 (□) and 2500 (■). (a) Ci(Z) vs ln(7) : some continuous cross-over is observed 
from a decorrelated regime at large scales down to a smooth one at small scales, (b) C2H) vs ln(Z) 
: a scaling behavior is obtained with the same slope — C2 = —0.025 ±0.003 for all signals, (c) C^{1) 
vs ln(Z) : in the inertial range and for the largest values of R\, the slope is compatible to zero up 
to finite sample effects, (d), (e) and (f) : same curves as in (a), (b) and (c), when all integral scales 
are set to 1. In the inset of (d), the solid line corresponds to the prediction of the Langevin model 
defined in the text after adjusting the integral scale L and the parameter a to the measured ones. 
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FIG. 2: (a) Magnitude connected correlation functions Cu(Ax) =< ln(\6vi(x)\) ln(|<5^(x+Ax)|) > c 
vs ln(Ax). (b) Square root of Cu(Ax) vs ln(Ax); Cu(Ax) behaves as a In 2 (Ax) whatever the 
experimental set-up with a systematic decrease of a when R\ is increased, (c) Same curves when 
all integral scales L are set to 1. The considered scale is I ~ 50 ^m for Castaing's signals, 0.15 mm 
for Baudet's signal and 3.2 mm for Gagne's signal. The symbols correspond respectively to the 
seven same signals as in Fig. 1. 
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FIG. 3: Prefactor a(R\) of the quadratic behavior of Cu(Ax) vs ln(Ax) (see Fig. 2). (a) ln(a) vs 
ln(Rx); the data are compatible with a power-law behavior a ~ R^ a with a = 0.20 ± 0.03 (solid 
line), (b) a as a function of l/ln(i?^); the data are compatible with a behavior a ~ K/\n{R\) 
with K = 0.027 ± 0.006 (solid line). The symbols have the following meaning : (o) Castaing's 
group, (x) Baudet's group and (■) Gagne's group. 
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